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Abstract: The primary physics signal events in the ICAL at INO are the y u charged current 
(CC) interactions with a well defined muon track. Apart from these events, ICAL can also detect 
other types of neutrino interactions, i.e. the electron neutrino charged current interactions and 
the neutral current events. It is possible to have a dataset containing mostly vyCC events, by 
imposing appropriate selection cuts on the events. The CC and the neutral current events form 
the background to these events. This study uses the Monte Carlo generated neutrino events, to 
design the necessary selection cuts to obtain a vyCC rich dataset. An optimized set of constraints 
are developed which balance the need for improving the purity of the sample and having a large 
enough event sample. Depending on the constraints used, one can obtain a neutrino data sample, 
with the purity of v e events varying between 55% to 70%. 
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1. Introduction 

India-based Neutrino Observatory or the INO, is an experimental facility to be set up in the southern 
part of India. An important component of INO will be the Iron Calorimeter (ICAL). The ICAL 
aims to study the interactions of atmospheric neutrinos and antineutrinos. The determination of 
the neutrino mass hierarchy is one of its prime objects, apart from adding to the precision of the 
oscillation parameters fHH - 

ICAL is a giant magnetized neutrino detector, with Resistive Plate Chambers (RPCs) as the 
active detector elements B.ELBD. Efficient tracking abilities, good resolution of energy and timing, 
good identification of the charge of the particles are the essential capabilities of this detector [0.0] . 
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ICAL comprises of 3 modules, with ~30,000 RPCs, and 151 iron layers weighing about 
50kton in total. Each module contains 8x8 RPCs in a layer, each of which spans a surface 
of 1.84 m x 1.84 m and is 2.5 cm in width (height) [0]. The steel structures to support them, how¬ 
ever occupy gaps of 16 cm in between the RPCs. An additional pair of slits is also present in each 
module, making way for the current carrying coils of the magnet. The RPC layers are interspaced 
with iron plates of 5.6 cm thickness, and an air gap of 4 cm, thus making the successive RPCs 9.6 
cm apart from each other in the vertical direction. 

The iron layers not only serve as the heavy target for most of the neutrinos but also cany 
the solenoidal magnetic field in the detector [0]. The neutrinos produce charged/neutral particles 
during interactions, which propagate through the detector. The iron layers are magnetised to ~ 1.3 
T. This accounts for the curvature in the tracks of the charged particles. The reconstruction of these 
tracks tells us the momentum and the charge of the particle. The kilometre-thick rock covering 
the INO cavern cuts down the cosmic ray particles. Hence, most of the events in the detector 
are produced by neutrino interactions only. Using the information of neutrino interactions in the 
detector, we can study the properties of these atmospheric neutrinos [0. 

The primary cosmic ray particles (like protons) interact with the atmospheric nuclei in the up¬ 
per layers. To a large extent, these interactions produce pions, which eventually decay into muons 
and muon neutrinos. The muons can further decay, leading to both muon neutrinos and electron 
neutrinos in the final state. Effects of oscillations can be seen on the upward going neutrinos which 
pass through the earth lUIUl . Oscillations lead to the creation of the tau neutrinos in the neutrino 
flux. 

The neutrinos and anti-neutrinos undergo two types of interactions, depending on the medi¬ 
ating particle. The charged current (CC) interactions are mediated by the W ± particles and the 
neutral current (NC) interactions by the Z°. The CC events are different for each flavour and can 
be distinguished by the charged lepton in the final state. If the charge of this lepton can be deter¬ 
mined, then a distinction between neutrino and anti-neutrino interactions can be made. The NC 
events have a neutrino/antineutrino in the final state hence it is not possible to distinguish neither 
the flavour nor the neutrino from anti-neutrino in these events El. 

The CC interactions can be quasi elastic (QE), resonance scattering (RS) or deep-inelastic 
scattering (DIS) type. The QE interaction gives a lepton corresponding to the neutrino flavor. The 
RS interaction produces one/two pions too in addition to the lepton. The DIS interaction, which 
dominates beyond the neutrino energy of 4 GeV, gives a shower of hadrons apart from the charged 
lepton. 

The NC interactions can be of RS or DIS type only. The RS and DIS interactions produce 
hadrons and a neutrino which goes undetected. So, only the hadron showers are visible to the 
detector in these cases. In the case of NC elastic scattering, because of no hadron shower in the 
final state, the event cannot be detected. 

The presence of the thick iron layers puts a lower threshold on the minimum energy of a de¬ 
tectable particle. Therefore, the detection probability of sub-GeV neutrinos is very less, due to the 
trigger criteria, which removes most of the random noise events. The neutrinos in the intermediate 
energy range of all the three flavors may interact with the detector. 

The ICAL is more efficient in the study of muon neutrinos in the GeV range. The presence of 
the magnetic field further enables us to tell apart the ji tracks from the /i 1 tracks in the detector. 
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However, apart from these well-recognizable muon track events, i.e., v^/v^CC interactions, the 
ICAL will also contain v e /v e CC interactions and also the NC interactions of all three flavors lOl. 
These interactions do not give clear muon tracks like the v u CC events. Therefore, apart from the 
muon track containing events, we should also focus on these “muonless” events, in order to extract 
maximum possible information derivable from the ICAL detector. 

The paper is organized as follows. In section 2, the details about the generated data which we 
have used in the study are explained, along with the details of how an event is detected and char¬ 
acterised. Section 3 deals with the application of various selection criteria which aim to increase 
vyCC purity in the sample. The criteria are explained in detail, along with their physics justifica¬ 
tion. Their effects on the total dataset are included in the corresponding subsections. Some of the 
criteria are logically applied on the total dataset to check if the percentage of NC events can be 
enhanced, in section 4. Section 5 discusses the background contribution of the v e CC and the NC 
events to the v u CC events sample. In section 6, we have studied the contribution of the muonless 
events in determining the neutrino mass hierarchy. Section 7 concludes with a summary of the 
results and a few important comments. 

Very few experiments have been able to study v e /v e events at higher energies. IMB, Kamioka, 
Super-Kamioka and very recently MINOS and T2K data are perhaps the only examples. ICAL is 
sensitive to neutrinos in the energy range of GeVs. Therefore, obtaining a source of data rich 
in electron neutrinos in this higher energy range is an important aspect for the whole neutrino 
community. One will be able to study the characteristics of these high energy vy/vys from the 
atmosphere. Other possible physics studies are contribution of these events in the determination 
of the neutrino Mass Hierarchy, and also utilization of the NC events in searches for the sterile 
neutrinos [0, mi l. 

2. Simulations 

The following analysis is done with neutrino events generated by Nuance neutrino generator. To 
reduce the influence of statistical fluctuations, we have generated 500 years of data. To begin with, 
we use normal hierarchy parameters. Simulations with inverted hierarchy parameters have also 
been checked, but not mentioned here, to avoid repetition. The generated events are then simulated 
in the ICAL detector using GEANT4 HT5H . 

2.1 Signal-Detection 

Charged particles produced by the neutrino interactions pass through one or more RPCs and gen¬ 
erate hits. These hits are our primary signals. The layer number of RPC gives the z-coordinate 
of the hit. The x and y-coordinates are given by the copper-strips of the pick-up panels which are 
orthogonally oriented at the top and the bottom of the RPCs. The number of strips in x-direction, 
with a signal, gives x-strips and similarly in y-direction gives y-strips. The maximum of the num¬ 
ber of x-strips or y-strips is defined to be the number of “strips-hit” in that layer. This number of 
strips-hit in a layer, when summed over all the layers which have received hits in that event, gives 
the number of strips-hit in that event. The hits distribution mentioned hereafter, refers to this value 
of the number of strips-hit in an event. 
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2.2 Types of events and their signatures 

Vertical and high energy muons travel through large number of layers before stopping/decaying. 
Therefore, the v^s which have high energy and are incident mostly along the vertical direction give 
hits in more number of layers, in case of CC interactions. In fact, the muons thus produced, form 
clear tracks in the detector and their momentum can be reconstructed. The curvature induced along 
the muon path due to the magnetic field, leads to the charge identification. On the contrary, the v^s 
which have lower energy or arc incident mostly along the horizontal direction 111611 are confined 
in less number of layers. They can hardly be distinguished from the hadron showers which also 
emanate at the event vertex. 

The v e CC events produce electrons which can give rise to em showers, but no track can be 
seen. The NC events have no charged lepton in the final state and hence have lower number of hits. 
They have only a set of hadrons in the final state and, in general, are indistinguishable from vyCC 
events. mm 

3. Selection Criteria to obtain v e CC rich Sample: 

We devise certain conditions to ensure that the selected event sample contains mostly vyCC events, 
with minimum possible background of y u CC and NC events. 

The simplest information available to us about an event is: (i) the number of hits that an event 
gives in the detector and, (ii) the number of layers it is passing through. Therefore, one first focuses 
on them. 

In order to understand the behaviour of the different neutrino events in the detector, we first 
have a look at the nature of hits distributions of all the three event types: the NC, the vyCC and the 
V U CC types, in different energy ranges of the incident neutrinos. 
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Figure 1. Hits distribution in the neutrino energy bins (from top to bottom in order): E v ={0.1,0.8} GeV; 
E v ={0.8,5.0} GeV; E v ={5.0,20.0} GeV; E v ={20.0,100.0} GeV, in case of the three types of neutrino events 
(from left to right in each row): Vj.CC; all NC (+V T CC); V' U CC event types. The x-axis gives the number of 
hits, and y-axis contains the counts of the events. 
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The distributions in figure 1 show that, for v u CC events, the number of hits is greatly enhanced 
with increasing energy. This increase is much less for v e CC events and hardly exists in case of the 
NC events. The figure also clearly suggests a lower threshold of ~10 hits to suppress a large 
fraction of NC events and low energy v e , y u CC events. With this cut, 12% of the total NC events 
is retained. The survival fraction for vyCC events is about 18%. So, out of the total set of survived 
data of the v e CC and the NC, more than 60% are v e CC events, as seen in Table [T|. The events 
containing high energy muons can be separated by restricting the number of layers or by being 
identified by the track reconstruction algorithm of ICAL. 

Below, we describe the selection criteria to generate event samples with maximum vyCC 
events. 


3.1 Hits and Layers 


The nature of the hits distributions is observed for all the three types of interactions mentioned 
above. A selection cut on the number of hits appears to be an effective way to obtain an events 
sample with a majority of vyCC events. The NC events peak at very low number of hits, unlike the 
CC events. So, it is difficult to select an event sample rich in NC events. 

The electrons/positrons, in general, travel a shorter distance than the hadrons. On the other 
extreme, the muons of the y u CC events travel through several layers. A primary observation of the 
layer distribution shows that a cut on the number of layers hit in an event is an effective criterion. 
Here the “layers” refer to the number of those layers which receive one or more hits in an event. 
However, the layer cut is a very sensitive cut, owing to the thick iron layer in between two RPC 
layers. 

As mentioned earlier, y u CC events either with low energy muons and/or in horizontal direc¬ 
tion, do not have identifiable muon track. Such events have been found to be a significant contrib¬ 
utor in the selected events sample. So, a separate count is maintained for them, in order to make it 
easy to derive a strategy to reduce them. The rest of the backgrounds, i.e., the vyNC, v^NC, v r CC 
and NC are all contained in the “others”. The v r CC events being pretty small in number, are not 
separately counted. 
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Figure 2. Distribution of number of layers for all events with more than 10 hits, for the 500years NH data 
in E v = {0.1,100} GeV. 

The v e CC events layer distribution, after a minimal hits cut (say 10), peaks around 5 while 
that of the V U CC events peaks around 10 layers. So, by selecting events which are confined in 5 
layers or less, we can reject v u CC events which are energetic and/or vertical. However, low energy 
V U CC events, especially those in the horizontal direction, do pass this cut and give rise to the events 
listed below. Various cuts on the number of hits and layers have been imposed on the set of events 
(£^={0.1,100} GeV). A few significant ones among them are tabulated in Table [T]. 


Selection Criteria 

v e CC 

others 

v^CC 

hits>0 

1106742 

1050814 

1682527 

hits >10 (only) 

202838 

130642 

672566 

hits >15 (only) 

97535 

69340 

445977 

hits >20 (only) 

52398 

42476 

314597 

hits>15; layers<4 

47711 

19390 

19875 

hits>15; layers<5 

68702 

32953 

36211 

hits>15; layers<7 

89614 

52550 

76194 

hits>10; layers<4 

125321 

56177 

62113 

hits>10; layers<5 

163807 

82717 

107350 


Table 1. Events counts after applying the selection cuts on the Geant output of the NH 500 years data in 
E v ={0.1,100} GeV. 


Very high number of hits are mostly due to v u CC events. So, an upper cut on hits can eliminate 
such events. 

Average hits per layer criteria: This is a criterion derived by combining the above two, but still 
may help in eliminating the track containing events. The muon tracks give mostly 2-3 hits in a 
layer. So, applying a lower cut on the average hits per layer (hpl) seems to be quite reasonable in 
rejecting most y u CC events. The hits per layer cut is useful in studying vertical, high energy muon 
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events. Here, given that the number of layers is constrained to be < 5, it is not very effective. The 
addition of this cut leads only to a marginal improvement. 

Applying the above and more similar cuts, it is seen that the fraction of vyCC events in the 
sample increases, but at the cost of sample size and fraction of vertical events. Hence, an optimized 
set of criteria needs to be chosen from the varied sets of cuts mentioned here. 

3.2 Distribution pattern of the hits in the layers: 

The cuts on the basis of hits and layers are indeed the simplest and very effective selection criteria. 
However, a number of various other parameters have also been studied, to ensure how much they 
can contribute to improving the purity of the vyCC events in the sample. (Please refer to INO 
internal note OT9I1 for details.) 

The behaviour of an vy-interaction is certainly different from the other two types of interac¬ 
tions, as far our physics knowledge is concerned. The presence of the electron/positron makes it 
stand apart from the NC interactions. The way electrons/positrons lose their energies in the de¬ 
tector is different from the way fi + //I do. The challenge is to utilize these characteristics in 
distinguishing vy events from y u events in the data from the ICAL detector. 

We have so far observed the overall picture of the number of hits and layers in an event. As 
shown earlier in this paper, these cuts have left us with enough hope to invest more interest in this 
study. So, it is required that we look deeper, to decipher those characterizable features of the vyCC, 
which can be readily realised from the ICAL data. 

3.2.1 Maximum Hits Difference 

The vyCC events contain em showers. These should have generated a huge number of hits, but 
most of them are absorbed by the thick iron layers. However, in some events the shower may start 
at the verge of the iron layer. In such cases, a sudden and significant increment in the number of 
hits in the following layer is expected. 

The difference in the number of hits in two adjacent layers in an event is calculated. This 
difference is maximized over all such pairs in that event. The value of the maximum difference in 
hits thus obtained forms our present selection criterion. The effect of this cut is shown in Table 0. 


Selection Criteria 

VeCC 

others 

VuCC 

h>10; L<5; 

163807 

82717 

107350 

h>10; L<5; max hits diff.>5 

82500 

34701 

38824 

h>15; L<5; 

68702 

32953 

36211 

h>15; L<5; max hits diff.>5 

50295 

21844 

23991 


Table 2. Events counts after applying the hits-layers selection criteria and adding the cut on maximum 
difference in the number of hits in adjacent layers. (500 years NH data in E v ={0.1,100}GeV.)[“h”=#hits; 
“L”=#Layers ] 

This selection criterion adds to improving the vyCC events ratio by about 3-4%. However, 
a simultaneous study of the known (Nuance) information shows that the larger hits difference is 
given by mostly horizontal vyCC events. 










3.2.2 Comparing the hits in each layer: 


The number of hits in every individual layer in an event is studied. This criterion, in a way seeks 
a pattern in the number of hits in adjacent layers. A variety of patterns are assumed and checked 
with the set of events. Two of them are stated below. The underlying logic still rests on the concept 
of the em shower. 

Additional hits in the next layer: One of the layers hit is chosen and additional 5 or 6 hits 
are demanded in the very next layer. All the layers in that event are also checked. The event to be 
selected must have at least one such a pair of layers. It is to be simultaneously noted that a lower 
threshold of 2 layers becomes inherent. 

Majority of hits in one layer: One can call this criterion a modified version of the earlier 
one. According to this criterion, the event must contain 50% or 60% of the total number of hits in 
a single layer. Therefore, no lower cut on the number of layers is required here. 

The effect of the selection cuts are tabulated in Table 0. 


Selection Criteria 

v e CC 

others 

VuCC 

v e CC purity % 

hits >15; layers <5; 

68702 

32953 

36211 

50 

hits>15; layers<5; h L >h L ±i+5 

47009 

21191 

22934 

52 

hits>15; layers<5; h L >50% 
hits 

38479 

13745 

16934 

56 

hits>15; layers<5; h L >60% 
hits 

29123 

9038 

11948 

58 

hits>15; layers<4; 

47711 

19390 

19875 

55 

hits>15; layers<4; h L >h L ±i+5 

34399 

13308 

13868 

56 

hits>15; layers<4; hz, >50% 
hits 

32737 

10931 

12679 

58 

hits>15; layers<4; h t >60% 

hits 

26006 

7735 

9690 

60 


Table 3. Events counts after applying the hits-layers selection criterion and demanding (i) 5 additional 
hits in adjacent layers (h/ , hz.±i); (ii) 50-60% of total number of hits in one layer. (500 years NH data in 
E v ={0.1,100} GeV.) [“hits”=total #hits; “h/ ”=hits in any of the layers, say the V 1 ' layer.] 


3.2.3 The Overall Distribution Pattern of Hits from the average number of hits among the 
layers 


The hits in different layers of the v e CC events are non-uniform. The hits are mostly over concen¬ 
trated in some layers, while entirely sparse in the rest. 
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Figure 3. Schematic Diagram of hits in the RPC layers 

Figure 3 shows, the hit pattern among various layers in an event (left panel) and the number of 
hits vs layer number (right panel). For the plot in the right panel, the lowest layer hit is labelled to 
be 0, the next layer is 1 and so on. In such a plot, the v u CC gives a broader peak than the v ( ,CC / 
NC. Hence, events selected with such sharper peaks should reduce the fraction of v u CC events in 
our sample. This property can be parametrized as either the mean or RMS value of the layerwise 
hits distribution of each event. However, having studied both the quantities, the cut on the RMS 
value appeared comparatively more effective. Some of the results are shown in Table [|. 


Selection Criteria 

v e CC 

others 

v M CC 

v e CC pu¬ 
rity % 

h>15; L<5; 

68702 

32953 

36211 

50 

h>15; L<5; rms<1.2 

56254 

24916 

25431 

53 

h>15; L<5; rms<1.2 ; max hits diff.>4 

48248 

20452 

21241 

54 

h>15; L<5; rms<1.2 ; max hits diff.>6 

39610 

15585 

16969 

55 

h>10; L<4; 

125321 

56177 

62113 

51 

h>10; L<4; rms<1.2 

111858 

47961 

52860 

53 

h>10; L<4; rms<1.2; max hits diff.>3 

86157 

35115 

37026 

54 

h>10; L<5; rms<1.2; max hits diff.>3 

99814 

43409 

46455 

56 

h>10; mean<2; rms<1.2; max hits 
diff. >3 

83954 

35130 

36127 

54 

h>10; mean<2; rms<1.2; max hits 
diff. >5 

60959 

23063 

24129 

56 

h>10; mean<2; rms<1.2; max hits 
diff.>5; hpl>4 

51249 

18247 

18922 

58 


Table 4. Events counts after applying the hits-layers selection criterion; adding the cut on the variance from 
the mean of the vertical distribution of hits in layers, i.e. rms; the criteria of max hits diff. included for 
further improvement. (500 years NH data in E v ={0.1,100} GeV.) [“h”=#hits; “L”=#Layers; “hpl”=avg 
hits/layer.] 

The RMS cuts appear quite effective in improving the ratio. In fact, the cut of maximum 
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difference in the hits further adds to bettering the results. Therefore, one can obtain ~55% majority 
of v g CC events, with a moderately large sample size. 


4. The NC events fraction 

Our primary/main focus in this paper is to obtain a v e CC rich sample of neutrino events at ICAL. 
So, the selection cuts so far have been favoring vyCC events. However, the NC events fraction can 
also be enhanced comparatively. 

The cuts are based on the simplest criteria of hits and layers. However, the extent of this 
possibility depends on the hardware threshold to be put in the ICAL, for accepting an NC event. 
This requires a trigger algorithm different from that for the muon track-containing events. 


xIO 3 



#hits 

Figure 4. Distribution of number of hits for all non-zero hit events with E v = {0.1,100} Ge V for the 500years 
NH data. 
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Figure 5. Distribution of number of layers which received one or more hits in an event, E v = {0.1,100} GeV 
for the 500years NH data. 
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The NC events give mostly very less number of hits, and are confined in very few layers, as 
in figure 4 and 5. So, one might be tempted to put an upper threshold on the number of hits, to 
obtain an NC events rich sample. But, to deal with such a small number of hits in an event, one 
must be aware of neutrinos with sub-GeV energies too. So, the dataset including neutrinos with E v 
= {0.1,100} GeV is satifactory to proceed further. 


Selection Criteria 

v e CC (NH 

evts) 

others (NH 

evts) 

v^CC (NH 

evts) 

NC 

purity % 

0<hits<10 (only) 

903904 

920172 

1009961 

32 

0<hits<10; layers<2 

659926 

724065 

478480 

39 

0<hits<4 (only) 

406705 

568177 

345322 

43 

0<hits<4; layers<2 

397895 

558109 

321648 

44 

0<hits<4; layers = 1 

287799 

436263 

198539 

47 

4<hits<10; layers = 

1 

70330 

38360 

31046 

27 


Table 5. Enhancing NC fraction: Events counts after applying the selection cuts on the Geant output of the 
NH 500years data files in E v ={0.1,100} GeV 


The NC counts in the selected sample are almost equal to the sum of the selected vyCC and 
the v u CC events. The NC events have a very small number of hits in general. The cuts used earlier 
demanded a minimum of 10 hits and hence discriminated against the NC events. To favour NC 
events, one should design a cut demanding a very small number of hits. The dominance of the NC 
events can be gradually realized in case of events having 10 hits or lower, as in Table 0. 

If we demand # hits < 3 and # layers = 1 or 2, we get the event samples shown in the above 
table, which are quite rich in NC events. If the noise is kept under control, such events can be used 
to study mixing with sterile neutrinos. Trigger efficiency will play a major role in selecting such 
events. In fact, it has been checked that a sample of single hit events has more than 50% NC events 
in the sample. But obviously, just one-hit is an unacceptable criteria. Therefore, the selection cuts 
will have to be redesigned entirely, to obtain an events sample containing NC events in significant 
majority. 

5. v e CC and NC events as background to v u CC events 

Muon track containing events are the primary data for the ICAL, especially those within the range 
E v ={0.8,20}GeV. The v^CC events detected at ICAL must pass through (i.e. give hits in) a min¬ 
imum number of layers (5 or 6), so that the muon track can be reconstructed. This layer cut will 
undoubtedly select mostly the v u CC events, which form the signal events in this case. However, 
some v e CC and NC events also will pass through this cut and form background to this events 
sample. 

As shown in Table |6j, out of all the generated events, about 20% of the “others” do not give 
any hits in the ICAL. For v e CC and y u CC events, this fraction of “undetectable” events is about 
5%. The layers distribution of these events (for the energy range E v ={0.1,100}GeV) is shown in 
figure 5. 
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Selection Criteria 

V e CC 

others 

(NC+v t CC) 

Vju CC 

all generated events 

676014 

820854 

1103263 

#events with hits>0 

in ICAL 

649487 

678590 

1087709 


Table 6. Events counts before applying the selection cuts on the Geant output of the NH 500years data files 
in E v ={0.8,20}GeV. 

Since reconstructable v^CC events demand a minimum number of layers to be hit, the distri¬ 
butions with two such layer-cuts are quantified in Table [7j. These basically fetch us the number of 
events or percentage composition of the selected events sample. This feature of large suppression 
of the v e CC and NC events with this cut is evident in figure 5. 


Selection Criteria 

V e CC 

others (NC+v T CC) 

CC 

#events: L>5 

84115 

73849 

683635 


~10% 

~ 9% 

OO 

2 

#events: L>6 

35678 

37031 

579760 


~5% 

~6% 

~89% 


Table 7. Events counts after applying the selection cuts on the Geant output of the NH 500years data files 
in E v ={0.8,20}GeV. 


6. Contribution of the muonless events to v Mass Hierarchy 

In this section, we study effect of muonless events in the mass hierarchy determination. This is 
expected to be much smaller compared to that on the muon events. Nevertheless, we pursued this 
study with the hope of improving the hierarchy sensitivity of ICAL. 


6.1 Physics Motivation and Application: 

The matter effect modifies neutrino oscillation probabilities. For long pathlengths (L > 5000 km) 
and moderately large energies (5 GeV < E v < 10 GeV), matter effects lead to large changes in 
both P(Vp —> Vfi) (Puu) and P(v e —t v^) (P e n)- These changes can lead to an observable change 
in the muon event rate. By measuring this change, it is possible to determine the neutrino mass 
hierarchy. The oscillation probabilities involving v e , P(v e —>• v e ) (P ee ) and P(v e ) (P^ e ), also 
undergo large changes due to matter effects. The spectrum of the electron events is given by 


dN e 

dE v 


d&e 

dE v 


Pee + 


d<t> f i 

dE v 





Since muon neutrino flux d<t>u/dE v is twice the electron neutrino flux d<X> ( ,/dE v and the change 
in Pi le is half the change in P ee , the effect of these large changes mostly cancel each other out in 
the electron event sample. This fact makes finding matter effects in muonless events even more 
challenging. 
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6.2 The Generated Events Sample: 

The data files from Nuance, in the energy range G v ={ 0.1.100 }GeV arc fed into the Geant4 INO 
ICAL code to get the events sample for the following studies. The neutrino oscillations have been 
applied using the normal and the inverted mass hierarchy parameters, which are denoted as NH 
and IH respectively. The oscillation parameters used are as follows: A/M 21 2 = 7.5 x 10 5 eV 2 , 
\Am e ff 2 1 = 2.47 x 10~ 3 eV 2 , sin 2 812 = 0.31, sin 2 20i3 = 0.09, sin 2 823 = 0.5 and 8 cp = 0. 

6.3 The Average MH y 2 : 

Since we are using Nuance Monte-Carlo to calculate the number of events, we need to take into 
account the MC fluctuations. Thus, if we simulate the NH events twice, with two different seeds, 
the y 1 between these two event samples will be non-zero. In fact, such X[ ruc = X 2 (NH I — NH2) 
will be approximately twice the number of bins. In addition, we calculate ^f alse = y 2 (IH — NH). 
If the NH is the true hierarchy, then we expect ^ 2 alse to be aprreciably greater than Xtme- 

To minimize the overall effect of MC fluctuations, we do our calculations for very large statis¬ 
tics and scale them down to 10 years. Here we consider data for 500 years. We have simulated the 
data for NH with three different seeds and similarly for IH. Thus, we have six values of ^ 2 (true) 
and nine values of x 2 (false). We take the average of each and define the average y 2 for hierarchy 
as 

< Z 2 > = < Xfalse >-< xLe > ■ 

The numbers of the vyCC and the v u CC events generated for each of the seeds, are listed in Table |8|. 


Sample ID 

NH v e CC 

NH v M CC 

IH v e CC 

IH v^CC 

seed 1 

676014 

1103263 

671309 

1103667 

seed 2 

674971 

1103879 

670827 

1105891 

seed 3 

675963 

1102817 

669664 

1104746 


Table 8. Events counts of v e CC and v^CC out of the 500 years Nuance data files, before interacting with 
the ICAL detector. Here, only the energy range £V={0.8,20} GeV is mentioned, which makes the major 
contribution in the value of y 2 . 

6.4 Calculation of Average y 2 assuming Normal Hierarchy (NH): 

The events are simulated in the energy range E v ={ 0.1,100} GeV, for both NH and IH, each with 
three different seeds. We use the earlier mentioned criteria to select tyCC rich samples in each 
case. To compare the distributions of these events in different cases, we sort them into a number of 
bins. We consider four different binning schemes. They are: 

• 1-bin scheme: The events are all contained in one single bin. Each of these events must have 
a minimum of 11 hits and are confined in 4 or less number of layers. 

• 3-bin scheme: The selected sample is divided into 3 bins, based on the number of hits. The 
events in the first bin should have a minimum of 11 hits but < 20 hits. The second bin 
covers the range of 21 to 40 hits, while events with 40 to 100 hits are put in the third bin. 
Appropriate layers or hpl cuts are also included in each of them, so that the sample contains 
maximum number of v e CC events. 
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• 10-bin scheme: The events are classified into 10 uniformly divided bins in the hits range 11 
to 100, along with suitable layers or hpl cuts. 

• 15-bin scheme: The hits range 11 to 20 is divided into 10 bins. The events giving hits from 
21 to 40 are grouped under 4 uniform bins. The fifteenth bin comprises of the events giving 
more than 40 hits. The layers or hpl cuts are imposed to ensure the majority of the v e CC 
events in the sample. 

We expect the < % 2 > to increase with the increasing number of bins, until a saturation value 
is reached. 

The values for x} r ue are calculated and tabulated in Table |9|, where as Table contains the 
values of xj aIse . 


Sample pairs 

xr (i) 

xr (3 ) 

* f 2 (10) 

* 2 (15) 

NH2-NH1 

2 

6 

10 

30 

NH3-NH1 

1 

3 

9 

27 

NH3-NH2 

1 

9 

18 

17 

IH2-IH1 

3 

9 

14 

44 

IH3-IH2 

0 

14 

23 

22 

IH3-IH1 

2 

7 

19 

47 

Average 

2 

8 

16 

31 


Table 9. Values of Xu-ue or X? f rom the three possible combinations of the three sets of dataset assuming 
Normal Mass Hierarchy for the binning schemes: (a) 1 bin-scheme; (b) 3 bin-scheme; (c) 10 bin-scheme; 
(d) 15 bin-scheme. 


Sample pairs 

Xf (i) 

X?(3) 

x}( 10 ) 

X 2 f ( 15 ) 

NH1-IH1 

44 

55 

72 

82 

NH1-IH2 

26 

29 

46 

62 

NH1-IH3 

39 

39 

50 

70 

NH2-IH1 

27 

29 

43 

50 

NH2-IH2 

13 

14 

31 

63 

NH2-IH3 

23 

28 

39 

54 

NH3-IH1 

35 

53 

80 

74 

NH3-IH2 

19 

23 

45 

67 

NH3-IH3 

30 

33 

50 

60 

Average (all) 

28 

34 

51 

65 


Table 10. Values of X/aise or Xf f r °m the nine possible combinations of the three sets of dataset assum¬ 
ing/expecting Normal Mass Hierarchy, (a) 1 bin-scheme; (b) 3 bin-scheme; (c) 10 bin-scheme; (d) 15 
bin-scheme. 


In Table 


11 


we have listed < x 2 >=< Xfalse > ~ < Xu-ue > an d the standard deviation in 
< r 2 >• This standard deviation is simply the sum of the standard deviations from the mean values 
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of xlue and Xf alse - 


hits binning 

< X 2 > 500yrs 

a(< x 2 >) or ^|500 yrs 

< X 2 > in 10 yrs 

1 

26 

10 

0.5 

3 

26 

16 

0.5 

10 

35 

20 

0.7 

15 

34 

20 

0.7 


Table 11. Values of average % 2 with standard deviation as the error in varying number of bins. The values 
are also scaled down to 10 years. 


A cos 0 binning has been attempted, but it appears to be of no additional help in the present 
study. The binning in terms of number of hits seems to suffice. 

In the above calculations, we assumed that NH is true. The results for the case IH is true are 
very similar. 


6.5 The Frequentist Approach with 10 years events samples 


The statistical fluctuations in the data are intrinsic to it. The Monte Carlo simulations mimic these 
fluctuations but every different simulation carries with it a different set of fluctuations. To properly 
estimate the effect of statistical fluctuations on < x 2 >> one should do a number of NH and IH 
simulations. From these, one can obtain the average value and the standard deviation of Xtrue an d 
X false • The difference in the average values of Xtrue ar| d X false c' ves the expected mean value of the 

< Z 2 > and the sum of their standard deviations gives us the probable range of < x 2 > EDD- 

We have divided the 500 years of data into fifty 10 years data samples, for both NH and IH. 
Using these, we computed 1225 values of NH-NH Xtrue ar| d 2500 values of NH-IH Xfalse- Their 
distributions are given in figure 6. The < xj a he > ' s higher than the < Xtrue > hy almost 3 (i.e. 

< X 2 > ~3). 
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Figure 6. Distribution of ^f alsc & Xu-ue calculated with lOyears of NH data and IH data, assuming NH 
ordering to be true. 


The plot in figure 6 show the Xmse ar *d Xtme distributions in the 15-bin scheme. However, both 
the Xfaise an d the Xtme have very broad distributions. Hence, one is not able to make any conclusive 
statement on the < % 2 >. 


7. Summary and Conclusion 


The hits and layers cuts have been confirmed to be the most important criteria for selecting a v e CC 
rich events sample. They alone fetch us an event sample containing around ~ 50% v,,CC events. 
The selection criteria to be finally chosen depend on the requirements of the physics study. One 
might insist on the maximum possible purity of the V,,CC events, even compromising the vertical 
events fraction or small sample size. 

The effects of the various other additional selection criteria are tabulated in Table P . 
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Selection Criteria 

Best Ratio of 

#v e CC to total 

#events 

Max. sample 
size with the 
ratio (k=10 3 ) 

Remarks 

Maximum Hits diff. 

53% 

156k 

Large sample size 

Comparison: hits in layers 

60% 

43k 


Overall Pattern: hits in lay¬ 
ers 

62%(58%) 

26k(88k) 


Single layer hits 

68% 

6.5k 

Very small sample 

size; Single layer 

more prone to noise 


Table 12. Effects of different selection criteria on the event sample after putting cuts on #hits and layers. 
Purity Ratio of v e CC in the total sample, and the corresponding sample size for 500 years NH data are 
shown. 

The results of optimizing the selection criteria may also be summed up in the following man¬ 
ner: 


• Maximum obtainable v,,CC purity: ~ 60% with counts ~ 100 events per year or ~ 62% 
with counts ~ 50 events per year. 

• Maximum obtainable vNC purity: ~ 47% with counts ~ 1800 events per year, provided 
there is no noise. 

• The selection criteria described in this report retain a majority of horizontal or near horizontal 
events. Any efforts to increase the percentage of vertical or near vertical events requires a 
compromise on the paid of purity of the v e CC events. 

A Nuance based analysis of the selected events has been done to understand the composition 
of the selected events. Three types of interactions are studied. The v e CC events chosen by selection 
cuts are thrice the background in the QE interaction bin. The NC and the v u CC sum up to just half 
the number of v e CC events in the RS bin. The DIS events count of the v,,CC almost equals the sum 
of other two. This implies, that the selection cuts will distinguish the QE and RS v e CC events over 
the rest, but not the DIS events. But the number of QE events are very less in comparison to the 
other two types. So, the overall signal to background ratio is not so good as can be expected. 

In the selection of y u CC events sample, V,,CC + NC events form a 20% background, if the 
events have a minimum of 5 layers hit. If this minimum is raised to 6 layers, the background 
contribution comes down to 10%. We also conclude that the contribution of the muonless events in 
determining the neutrino mass hierarchy is not zero. But the statistical fluctuations in the data arc 
too large for this contribution to have a significant effect. 
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